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ABSTRACT 

We use the ROSAT All-Sky Survey to study the X-ray properties of a sample of 
625 groups and clusters of galaxies selected from the Sloan Digital Sky Survey. We 
stack clusters with similar velocity dispersions and investigate whether their average 
X-ray luminosities and surface brightness profiles vary with the radio activity level of 
their central galaxies. We find that at a given value of cr, clusters with a central radio 
AGN have more concentrated X-ray surface brightness profiles, larger central galaxy 
masses, and higher X-ray luminosities than clusters with radio-quiet central galaxies. 
The enhancement in X-ray luminosity is more than a factor of two, is detected with 
better than 6cr significance, and cannot be explained by X-ray emission from the 
radio AGN itself. This difference is largely due to a subpopulation of radio-quiet, 
high velocity dispersion clusters with low mass central galaxies. These clusters are 
underluminous at X-ray wavelengths when compared to otherwise similar clusters 
where the central galaxy is radio-loud, more massive, or both. 



1 INTRODUCTION 



There is increasing evidence that the majority of radio-loud 
AGN at low redshift may be triggered by the accretion of 
hot gas. U s ing a sample of 625 nearby groups and clusters, 
iBest et aU (|2007l ) showed that the galaxies located closest to 
the centres of the clusters are more likely to host a radio-loud 
AGN than other galaxies of similar stellar mass. lAllen et al.l 
l|2006l ) analyzed Chandra X-ray images of nine nearby X-ray 
luminous elliptical galaxies and showed that the jet power 
is tightly correlated with the Bondi accretion rate onto the 
central black hole as estimated from the observed gas tem- 
perature and density profiles. 

It is also now recognized that the 'over-cooling prob- 
lem' at the centre of many galaxy clusters could be solved if 
radio-loud AGN can heat surrounding gas. Direct evidence 
of AGN heating came with the discovery of X-ray cavi- 
ties in the hot intracluster medium of a number of clusters 



and groups (iBoehrinEer et al.l Il993l: 


Churazov et all l2000l: 


Fabian et al.l 2000l; McNamara et al. 


2000l: Blanton et all 


20031: iGitti et al.ll2007h. Some of these cavities coincide with 



extended lobes of radio emission produced by an AGN in 
the central cluster galaxj0. The observations suggest that 
relativistic radio plasma displaces the X-ray emitting gas. 



generating turbulence and wave activity which heat the in- 
tracluster medium (ICM). 

AGN feedback may also explain why low-temperature 
clusters are less X-ray luminous than predicted by a homolo- 
gous scali ng of the properties o f hotter and more massive sys- 
tems (Nath fc Rovchowdhurvi,2002. : Best et all |2007| ') . The 
energy available from the central AGN depends primarily on 
black hole accretion rate. In low-mass groups, it can be com- 
parable to the total gravitational binding energy of the X-ray 
emitting gas, whereas high-mass clusters will be less strongly 
affected. As a result, relations such as Lx —Tx (X-ray lumi- 
nosity vs temperature) and Lx — o" (X-ray luminosity vs ve- 
locity dispersion) may be modified by radio-source heating. 
At smaller scales, AGN feedback may play an important role 
in the formation and evolution of galaxies. Recent work has 
demonstrated that it can plausibly explain the exponential 
cutoff at the bright end of the galaxy luminosity function, 
as well as the 'down-sizing' of galactic star-formation activ- 
ity at recent cosmic epo chs ()Croton et al. 2006; Bowo rTt al.l 
l2006l : iKang et al.|[2006l ) . 

It may be that hot gas accretion, AGN triggering, re- 
heating of ambient gas, and suppression of AGN activity 
occur in a cycle. Accretion onto the central black hole may 



^ Some cavities do not h ave associated radi o emission above cur - 
rent detection limits [e.g. lEttori et al.l l l2002h , IClarke et al.l ||2005|') 



and ljetha et al~ 1 1 |200^ 1 but these hot bubbles may have been : 
flated by a previous cycle of nuclear activity. 



© 0000 RAS 



2 Shen et al. 



cause an outburst which removes the fuel supply for further 
radio activity. Without AGN feedba ck, the ICM reverts t o 
the state that triggered the outburst (jChurazov et al.ll2005l l. 
Although this picture is attractive, it remains to be ver- 
ified in detail. We do not know the necessary conditions 
to trigger or to quench radio activity. We also do not un- 
derstand the extent to which the global X-ra y properties 
of the gas are modified by AGN feedback(e.g. Rizza et al.l 
I2OOOI : iQmma fc BinnevI 12004 iNusser fc 811^120081 ). Onlv a 
few nearby clusters have deep enough X-ray images to re- 
veal low-density cavities and permit the energetics of the 
gas to be studied in a spatially resolved fashion. If we wish 
to study how the global state of the X-ray gas is linked to 
AGN activity in the central galaxy, we are forced to adopt 
a more statistical approach. 

It has long been known that the radio activity of clus- 
ter central galaxies is linked with the cooling flow or "cool 
core phenomenon l|Burnij|l99(il : iFabianlllOM ). This is the 
fact that many but not all massive clusters show a strong 
central peak in X-ray surface brightness, almost always co- 
incident with the brightest cluster galaxy (BCG), within 
which the directly inferred cooling time is less than the 
Hubble time. The fraction of cool-core clusters with a radio- 
loud BCG is much higher than the fraction in the rest of 
the population. Furthermore, because the cores add sig- 
nificantly to cluster luminosity and reduce the emission- 
weighted temperature, cool-core clusters fall systematically 
to one side of many of the X-ray scaling r elations for clusters , 
e.g. the Lx — Tx and Lx — o" relatio ns (jFabian et al.l 1 19941 : 
lO'Hara et all I2OO6I : IChen et all 120071 ). Thus, in clusters of 
given velocity dispersion (or mass), cool cores tend to be as- 
sociated with enhanced X-ray luminosity, with massive cen- 
tral BCGs, and with radio-loud BCGs. Given that the prob- 
abilit y of radio activity increases with BCG mass CBest et al.l 
I2OO5I ). that both BCG mass and cluster X-r a,y lum inos- 
ity increase with cluster mass (|Edge fc StewartHl99ll ), and 
that central black hole ma ss correlates tightly with galaxy 
mass l|Tremaine et al.ll2002l ). it seems clear that radio activ- 
ity is tightly related to both black hole mass and the pres- 
ence of a cooling hot atmosphere which can provide fuel. 
The exact causal relation between these phenomena remains 
unclear, however, and recent observational studies of lower 
mass systems appear, at least superficially, in conflict with 
the tren ds established for re lative ly massive c l uster s [e.g. 
compare ICroston et~all (|2005l ) and Ijetha et akl l|2007h with 
the above references]. 

In this paper, we provide improved statistics for the 
relationship between cluster properties, both X-ray and op- 
tical, and the radio activity of their BCGs. We use a low 
redshift [z < 0.1) sample of 625 groups and clusters with 
ca refully controlled BCG identifications. These were selected 
bv lvon der Linden et al. (2007) from the Sloan Digital Sky 
Survey(SDSS, York et al.il200d ). Radio prope r ties of these 
BCGs were determined following iBest et al] l|2005l ), who 
cross-matched the galaxies with the National Radio As- 
tronomy Observator y Very Large Array Sky Survey (NVSS, 
ICondon et al.l Il998l ) and the Faint Images of the Radio 
Sky at Twenty Centimeters Survey (FIRST, iBecker et al.l 
1 19951 ). By combining these two radio surveys, the identi- 
fication of radio galaxies is both reasonably complete (~ 
95 per cent) and highly reliable (~ 99 per cent). Radio- 
loud AGN are then separated from star-forming galaxies 



with det ectable radio em ission on the basis of 4000A break 
strength l|Best et al.ir2005l ). 

The sky coverage of this cluster sample is ~ 4000 
square degrees. The only X-ray survey with enough sky 
coverage to provide a reasonable match is the ROSAT 
All-Sky survey (RASS). Because the RASS is quite shal- 
low, relatively few individual groups and clusters have suf- 
ficient X-ray flux for unambiguous detection, particularly 
at low velocity dispersion. The mean flux of such objects 
can be detected, however, by stacking their X-ray images 
(e.gjBartelmann fc Whit3l2003l : IShen et al.|[2006l : iDai et all 
I2OO7I : iRvkoff et al.ll2008l ). Such stacking avoids selection bi- 
ases that can occur if one analyses samples in which only 
the most X-ray luminous systems are detected. The limited 
resolution of RASS maps and our stacking strategy, make it 
difficult to study the structure of the ICM in detail, but we 
will see that concentration differences, refiecting the pres- 
ence of cool cores, are nevertheless detectable. 

Our paper is organized as follows. In Section 2, we in- 
troduce our sample of groups and clusters and describe their 
radio properties. In Section 3, we describe our X-ray detec- 
tion techniques, both our method for detecting individual 
groups and clusters and our stacking technique. In section 
4, we study and compare the X-ray properties of clusters, 
emphasising how the Lx — o" relation and the surface bright- 
ness profiles of clusters vary with the radio properties of their 
BCGs. We discuss the contribution of the radio AGN to the 
total X-ray emission in Section 5.1, while in Section 5.2, we 
discuss how the Lx — o" relation of clusters depends on the 
stellar properties of their BCGs. We present our conclusions 
in Section 6. 



2 SAMPLE 

We use the sample of groups an d clust ers of galaxies de- 
scribed in Ivon der Linden et al.l l|2007h , which is drawn 
fr om the C4 cluster catalogue of the SDSS Data Release 

3 (|Miller et al.ll2005l ). The clusters lie in the redshift range 
0.02 ^ 2 0.1. Von der Linden et al. developed improved al- 
gorithms for identifying the brightest galaxy (the BCG) and 
for measuring the velocity dispersion in each group or clus- 
ter. The velocity dispersion algorithm is designed to limit 
the effects of neighbouring groups and clusters. Clusters and 
groups with very few galaxies are also discarded. The veloc- 
ity dispersion a is measured within the virial radius -R200- 
The rad io properties of th e BCGs are taken from the cata- 
logue of iBest et all l|2005h . which has been updated to the 
SDSS data release 4. 

The final sample includes 625 groups and clusters, with 
velocity dispersion spanning the range from ~ 200 kms~^ 
to over 1000 kms~^ (see von der Linden et al. 2007 for more 
details). 134 out of the 625 BCGs have radio fiuxes larger 
than 5 mjy and are identified as radio-loud AGN. Five BCGs 
have radio fluxes which exceed 5 mJy but are clearly a result 
of star formation activity. There are 433 BCGs without a 
radio source brighter than 5mJy. The radio properties of 
the remainder are unknown because they lie outside the area 
covered by FIRST. Henceforth, we will refer to a group or 
cluster as 'radio-loud' if its BCG has been identifled as a 
radio-loud AGN, and 'radio quiet' if it is known to contain 
no radio source brighter than 5mJy. There are 134 radio-loud 
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Figure 1. Velocity dispersion and redshift distributions for all 
radio-loud clusters (solid histograms), for all radio-quiet clusters 
(dotted histograms) and for our matched control sample of radio- 
quiet clusters (dashed histograms). 



catalogues are not optimized for extended sources. An ex- 
tended source is likely to be deblended into several pieces 
by the ML search algorithm. 

There are carefully selected samples of X-ray clusters 
constructed from the RASS data, e.g. the Northern ROSAT 
All-S ky Galaxy Cluster Survey (NORAS, iBohringer et all 
|2000!) and the ROSAT-E SO Flux-Limited X-R ay Galaxy 
Cluster Survey (REFLEX. fBohringer et"aLll2004l )'). The RE- 
FLEX clusters, which have a completeness of ^ 90% at the 
flux limit of 3 X 10~^'^erg cm~'^ s~^, are mainly located in the 
southern sky {S ^ 2.5) and have little overlap with our SDSS 
clusters. For the NORAS clusters, the completeness is esti- 
mated to be ~ 50% at an X-ray flux of 3 x 10~^'^erg cm~^ s~^ 
l|Bohringer et aLll2004) . 

Our sample of groups and clusters has been selected 
from an optical galaxy catalogue. The centre of the cluster 
(taken to be the position of the BCG) and the virial ra- 
dius of each cluster have already been determined. We use 
this extra information when studying their X-ray proper- 
ties. We first check whether the position of the BCG of each 
cluster is consistent with a peak in X-ray emission (Section 
13. We then optimize th e growth curve analysis algorithm 
of Bohringer et all (|2004l ) to detect the X-ray properties of 
our clusters individually (Section l3.2|l . Finally, we develop a 
stacking algorithm to measure the mean X-ray properties of 
clusters as a function of their optical properties, independent 
of the detection limit of the RASS (Section 13. 3p . 



and 433 radio-quiet clusters in our sample. The fraction of 
radio-loud clusters in our sample is smaller than the fraction 
of clusters usually quoted as having cooling cores(~ 40% 
according to lPeres et alj[l998l ). This may well be a selection 
effect since X-ray selected cluster samples are biased towards 
massive and X-ray luminous systems. 

In figure [1] we show histograms of velocity dispersion a 
and redshift z for radio-loud and radio-quiet clusters. The 
two redshift distributions are very similar. The radio-loud 
clusters have slightly higher velocity dispersions than the 
radio-quiet objects. The median a of radio-loud clusters is 
428 kms~^, with 16 and 84 percentiles at 277 and 595 kms^^ 
respectively. For radio-quiet clusters, the corresponding a 
values are 392, 252 and 583 kms~^. The figure also shows 
a and z distributions for a control sample of radio-quiet 
clusters (dashed line) . This was constructed by choosing the 
radio-quiet cluster closest in a and z to each radio-loud clus- 
ter (see Section 14. 2p . This matching procedure was intro- 
duced in order to minimize possible observationally induced 
biases when comparing the two samples. 



3 X-RAY DETECTION 

We use the RASS to study the X-ray properties of our sam- 
ple of groups and clusters. The RASS mapped the sky in 
the soft X-ray band (0.1 — 2.4keV) with exposure time vary- 
ing between 400 and 40,000s, depending on direction on 
the sky. The resolution of the RASS images is 45 arcsec. 
Two source catalogues were generated based on RASS im- 
ages through a maximum-li kelihood (ML) sea rch algorithm: 
the bright source catalogue ( Vog es et al.lll99^ ) and the faint 
source catalogue (Voge s et a l. 2000l ). However, these source 



3.1 Finding the X-ray centre of the clusters 

Before the X-ray luminosity of a cluster can be measured, 
the position of its centre must be determined. In the optical, 
the position of the BCG is usually taken to define the cluster 
centre. However, this positio n may be offset from the peak 
of th e X-ray emission (e.g. iDai et all l2007l : iKoester et all 
l2007f) . In order understand whether such offsets are generic 
for clusters in our sample, we have developed an algorithm 
to determine the X-ray centre of a cluster. 

First, the ML search algorithm used to generate the 
ROSAT source ca talogue is applied to the RASS images 
(j Voges et alJI 19991 ). Sources with detection likelihood L > 7 
in the 0.5 — 2.0 keV band are retained. Although an extended 
source is frequently de-blended into several pieces by this al- 
gorithm, the central peak is always identified. We match the 
position of each BCG with the ML detections using a toler- 
ance of 5 arcmin in radius. X-ray sources found within this 
radius are considered as candidate X-ray centres. We find 
that 210 clusters have X-ray sources within a 5 arcmin ra- 
dius and some clusters have several candidate X-ray centres. 
We then exclude candidates which are clearly not associated 
with the clusters by eye. If a point source is detected within 
the 5 arcmin search radius but significantly offset from the 
BCG, we assume that the source is a contaminant and is 
not part of the cluster. We show such an example in Fig. 
Bl of Appendix B. For clusters that still have more than 
one candidate X-ray centre, we pick the closest one. Most 
of the sources we identified as the X-ray centres of clusters 
have extension likelihood greater than zero. This confirms 
that the X-ray emissions we identified are from the extended 
ICM rather than from point-like AGN. As we will show in 
Section [5. II the X-ray emission from radio AGN in the soft 
X-ray bands is expected to be much weaker than the X- 
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Figure 2. A histogram of the offset between BCG and X-ray 
centre. The upper panel shows the offset in units of arcmin, while 
the lower panel shows the offset in units of kpc. 



ray emission we measure. Our final sample of clusters with 
identified X-ray centres contains 157 objects. 

In Figure [2] we plot the offset between X-ray centre and 
optical centre (i.e. the BCG position) for these 157 clusters. 
The upper panel shows the offset in units of arcmin, while 
the lower panel shows the corresponding projected physical 
distance at the redshift of BCG. As can be seen, the X-ray 
centre is consistent with the BCG position in most cases. 
The typical offset is less than 2 armin 3 RASS pixels). 
This corresponds to a physical distance of around 50 kpc an d 
so is consistent with the results of Katavama et al.| (|2003l ) . 

Since the BCG positions are, in general, consistent with 
the X-ray centres, we will assume that the BCGs mark the 
centres of those clusters which are too faint to determine an 
X-ray centre directly with our algorithm. 



3.2 Measuring the X-ray luminosity 

After the centre of each cluster has been determined, we 
use a growth curve analysis to count the number of photon 
events as function of radius and so to evaluate the cluster 
count rate. 

The background of each cluster is estimated from a cen- 
tred annulus with inner radius -R200 and width 6 arcmin. 
To exclude contamination of this background estimate by 
discrete sources, we mask out all ML-detected [L > 7) 
sources inside this annulus. The surface brightness of the 
background is then calculated using the formula 



EN 
1=1 



(1) 



five area of the annulus after the contaminating sources are 
masked out. 

The count rates within -R200 may also be contaminated 
by discrete X-ray sources unassociated with intracluster X- 
ray emission. As above, we therefore mask out all the ML- 
detected (L > 7) point sources that are located more than 
2.25 arcmin (3 image pixels) from the X-ray centrfl As al- 
ready mentioned, an extended source is often detected multi- 
ple times by the ML algorithm. We therefore visually inspect 
all the ML sources and determine whether they are contam- 
inants or part of the X-ray emission from the cluster. We 
show an example of this process in Fig. B2 of Appendix B. 

The cumulative source count rate as a function of ra- 
dius is calculated by integrating the source counts in concen- 
tric rings and subtracting the background contribution. We 
integrate the source count rate to the X-ray extension ra- 
dius Rx- We determine this radius i n two different ways, as 
described in iBohringer (|2000l ). The first chooses the 

radius where the increase in source signal is less than the 
Icr uncertainty in the count rate. The other is the plateau- 
fitting method, where the slope of the plateau (in units of 
count rate per arcmin) is less than 1 percent. These two 
methods give consistent results for clusters where the count 
rate profile can be determined with high S/N. For low S/N 
clusters, we choose the better determination of Rx by vi- 
sually checking the cumulative count rate profile. A source 
is said to have a significant X-ray detection if the counts 
within Rx are three times larger than the Poisson fluctu- 
ation in the photon counts. We detect 159 galaxy clusters 
individually according to this criterion. 

To convert the count rates into X-ray fluxes, we assume 
that the X-ray emission has a thermal spectrum with tem- 
perature T, that the cluster gas has metal abundance equal 
to a third of the solar value (jRavmond fc Smithlll977l ) and 
that interstellar absorption c an be inferred from the G alac- 
tic hydrogen column density (iDickev fc Lockmanlll990l ). The 
gas temperature T is assumed to fo l low th e empirical scaling 
relation measured bv I White etal] ()l997[ ). 



f - V 

V403 kms-iy 



(2) 



where the sum is over the photon events, is the effective 
exposure time at each photon position, and S is the effec- 



One might question whether this assumption is robust. 
Most studies of the T—a relation flnd that it does not depart 
signiflcantly from the virial theorem ex pectation (T oc a^) , 
even for low-rnass galaxy groups (e.g. ICirardi et all Il996l : 
IWu et al.l 1 19991 : iMulchaevI I2OO0I : IXue fc Wull2000l ). We fur- 
ther note that a variation in temperature of 50 percent 
makes less than a 5 percent difference to the flux estimate 
for most of our clusters0 Throughout this paper we adopt a 



Here, we assume that all X-ray photons inside the 2.25 armin 
circle are emitted from the ICM. The probability that a fore- 
ground or background source is located inside this small circle by 
chance is close to zero. The X-ray emission from the central AGN 
itself is expected to be negligible compared with that from the 
ICM (see Section ISlTl. 

^ If the gas in radio-loud groups is hotter than in radio- quiet 
groups as suggested, for example, bv ICroston et al] fcOOSt) , we 
will underestimat the X-ray luminosity of radio-loud clusters and 
the conclusion we reach below about the difference in X-ray lumi- 
nosity between radio-loud and radio-quiet clusters (Section l4.2.2l l 
will be reinforced. 
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concordance ACDM cosmology with Hq = 70km s ^ Mpc ^, 
^0 = 0.3, and = 0.7. 

For some clusters, the total X-ray flux will exceed the 
flux measured within the extension radius Rx- To estimate 
the flux that is missed outside Rx, we adopt a /3- model 
surface brightness distribution Ix{R) with P — 2/3, 

/x(7^) = /o(l + U)-'^+°■^ (3) 

where Jo is the central surface brightness and Rc is the core 
radius and is assumed to be proportional to -R200, with Rc — 
0.14_R2oo for radio-loud clusters and Rc = O.I8-R200 for radio- 
quiet clusters The extension correction factor fs is then 
defined as 

C'"Ux{R)RdR 

fs = — — (4) 

J^'''' Ix(R)RdR 

The exact values of Rc and (3 are far from certain for each 
individual cluster. As a result, the X-ray luminosities of clus- 
ters with large correction factors fs have larger errors. To 
account for this effect, we assume that the correction factor 
/e has an uncertainty of 0.1/b. We note that for a few clus- 
ters with sufficiently high S/N, we can apply a /3 model fit 
and derive the extension correction factor /e- We find an- 
swers that are consistent with the estimates using Equ. Q 
with an uncertainty of ~ 10%. The error on our estimate of 
Lx for each cluster thus has two terms: the error of the flux 
estimation inside Rx and the error on the correction factor 

fE. 

The X-ray properties of the 159 clusters that are indi- 
vidually detected in RASS are listed in Appendix A. 

3.3 Stacking Analysis 

The clusters that are individually detected in the RASS are 
biased to the nearest and most X-ray luminous systems at 
each velocity dispersion. We can obtain an unbiased estimate 
of the average X-ray luminosity of our clusters by stack- 
ing objects with similar optical properties. We describe our 
stacking algorithm below. 

We divide the clusters into bins of velocity dispersion 
with a width of ~ 100 kms"'^. Within each velocity disper- 
sion bin, the angular size of the clusters varies substantially 
because of the spread in redshift. (Recall that -R200 is pro- 
portional to a.) We centre each cluster image on the X-ray 
centre if this has been identified (see Section [3. l|l and on the 
BCG otherwise. Before stacking the images, we scale them 
all to the same size in units of -R200 and mask contaminat- 
ing sources in the same way as in our analysis of individual 
clusters. The background is estimated as the average flux 
inside an annulus with inner radius -R200 and outer radius 
1.3i?200. The remaining steps in the detection algorithm (de- 
termining the extension radius Rx and count rate) parallel 
those used for individual clusters. A stacked image is said to 

^ Here, the different Rc values for radio-loud and radio-quiet clus- 
ters are based on the results from stacked images shown in Section 
14.2.11 This difference introduces a typical change of Lx less than 
5 percent. We would come to exactly the same conclusions if we 
adopted the same Rc relation for all clusters, regardless of their 
radio properties. 



have a robust X-ray detection if the number of source pho- 
tons within Rx is also three times larger than the Poisson 
fluctuation in the photon counts. 

For a stack of A'^ clusters with individual X-ray luminosi- 
ties Lx,i, redshifts Zi, average Galactic hydrogen columns 
Nh.z and average RASS effective exposure times ti, the num- 
ber of source photons that should be detected in the RASS, 
Ns, is 

N 

Ns='^ Lxj g{NH,i, Zi)ti , (5) 

where g{NH,i, Zi) is a function which converts the bolometric 
X-ray luminosity to observed count rate. Since the X-ray lu- 
minosities of the sources in a stack are similar, the weighted 
average X-ray luminosity of the stack can be defined as 

^x,s = . (d) 

L.i=i9{NH,i,Zi)U 

We fit the surface brightness profile of each stack with 
a P model to make the extension correction, and we correct 
the X-ray luminosity to the value expected within i?200 (see 
Section 1323. 



4 X-RAY PROPERTIES OF THE CLUSTERS 

In this section, we investigate how the X-ray properties of 
the clusters depend on the radio properties of their central 
BCGs. We focus on the comparison of the Lx — cr relations of 
radio-loud and radio-quiet clusters. We first show results for 
clusters that were detected individually in the X-ray images 
(Section 14. 1[) . We then present results for stacked cluster 
samples (Section 14. 2p . 

4.1 Results for individual X-ray detections 

4. LI Detected fraction 

As we showed in Section [3.21 only ~ 25 percent (159 out of 
625) of our clusters are individually detected in the RASS. 
The solid histogram in figure [3] shows the detected frac- 
tion as function of velocity dispersion. As expected, the de- 
tected fraction is higher for higher velocity dispersion clus- 
ters. It increases from less than 15 per cent for groups with 
a < 400 kms~^ to more than 50 percent for clusters in the 
highest velocity dispersion bin. 

We find that 55 of the 134 radio-loud clusters and 88 
of the 433 radio-quiet clusters are detected. The detected 
fractions as function of velocity dispersion for these two 
sub-samples are plotted in Fig. |3] As can be seen, the de- 
tected fraction is systematically higher for radio-loud than 
for radio-quiet clusters in all the velocity dispersion bins. 

4.L2 Lx — o" relation for clusters with individual X-ray 
detections 

In figure [4] we plot bolometric X-ray luminosity Lx as func- 
tion of velocity dispersion a for clusters with individual X- 
ray detections. Lx is in units of Lx/10''''erg cm~^ s~^(L44) 
and velocity dispersion in units of cr/500 kms~^((T5oo). 

We use the BCES orthogonal distance regression 
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Table 1. Fitting parameters for various Lx — o" relations, a and b are the fitting parameters of equation JT}. A'^ is the number of objects 
in each sample. For staclsed samples, A'^ is given as N1/N2, where Nl is the number of stacks while N2 is the total number of clusters 
in the sample. M is the median BCG stellar mass, b' gives the zero-points obtained if the Lx — f relations are refit requiring the slope 
a to be 2.97 for the low BCG mass and low BCG mass radio-quiet samples and to be 4.40 for the other stacked samples. 



Sample 






a 








b 




N 




M(M0) 


b' 


Individual, all 


4 


39 


± 





32 


-0 


530 ± 


037 


159 


1 


82 


X 


10" 




Individual, radio-loud 


4 


48 


± 





61 


-0 


535 ±0 


058 


55 


2 


09 


X 


10" 




Individual, radio-quiet 


4 


50 


± 





44 


-0 


530 ±0 


056 


88 


1 


78 


X 


10" 




Stack, radio-loud 


4 


40 


± 





53 


-0 


600 ±0 


099 


8/134 


1 


66 


X 


10" 


-0.600 ±0.032 


Stack, control radio-quiet 


4 


07 


± 





24 


-0 


935 ±0 


049 


8/134 


1 


35 


X 


10" 


-0.942 ± 0.045 


Stack, low BCG mass 


2 


88 


± 





17 


-1 


150 ±0 


055 


7/217 


7 


76 


X 


lOio 


-1.151 ± 0.033 


Stack, intermediate BCG mass 


4 


58 


± 





56 


-0 


712 ±0 


116 


8/208 


1 


35 


X 


10" 


-0.703 ± 0.037 


Stack, high BCG mass 


4 


16 


± 





44 


-0 


638 ±0 


072 


8/200 


2 


14 


X 


10" 


-0.646 ± 0.030 


Stack, radio-loud, low BCG mass 


4 


17 


± 





74 


-0 


572 ±0 


104 


6/67 


1 


20 


X 


10" 


-0.563 ± 0.041 


Stack, radio-loud, high BCG mass 


4 


49 


± 





47 


-0 


650 ±0 


072 


6/67 


2 


14 


X 


10" 


-0.651 ± 0.048 


Stack, radio-quiet, low BCG mass 


3 


07 


± 





21 


-1 


160 ±0 


053 


7/219 


8 


32 


X 


10" 


-1.162 ±0.032 


Stack, radio-quiet, high BCG mass 


4 


25 


± 





44 


-0 


704 ±0 


088 


8/214 


1 


74 


X 


10" 


-0.711 ±0.040 
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Figure 3. The fraction of clusters with individual X-ray detec- 
tions as a function of velocity dispersion. The solid histogram 
shows the result for the whole sample, while the dotted and 
dashed histograms show results for radio-loud and radio-quiet 
clusters, respectively. 



Figure 4. The Lx — o" relation of clusters that are individually 
detected in the RASS. The solid line shows the linear fit of equa- 
tion ((SJl. For comparison, the dotted line shows the fit obtained 
by Ortiz-Gil et al.{2004) for REFLEX clusters. 



method (jAkritas fc Bershadvl [l996l ) to fit a linear relation 
between log L44 and log crsoo , 



log L44 = a log (7500 + b. 



(7) 



This fitting method takes into account the observational 
errors on both variables and the intrinsic scatter in the 
relation. The determination of the error on Lx has been 
described in section I3.j1 w hile the error on a is given by 
Ivon der Linden et al. (|2007h . The best fitting relation is 



logL44 = (4.39 ±0.32) log (J500 + (-0.530 ± 0.037), 



(8) 



and is shown as solid line in Fig. |4] The fitting param- 
eters are also listed in Table [T] For comparison, we also 



plot the Lx — o" rela tion derived for REFLEX clusters by 
lOrtiz-Gil et"al] (|2004| ) using the orthogonal distance regres- 
sion method (dotted line). As we can see, the slopes of two 
fitting relations are consistent with each other within the 
1-cr uncertainty. The zero-point of our relation is slightly 
higher. This difference is caused by the fact that Ortiz-Gil 
et al. adopted a lower value of Rc for the ex tension correc- 
tion [Rc oc L^-^^, Equ. (6) of .Bohringer et al . (2000)]. 

We show Lx — ty relations for the sub-samples of detected 
radio-loud and radio-quiet clusters in the left and right pan- 
els of figure [S] respectively. The BCSE orthogonal regression 
method is again used to fit a linear relation between log L44 
and log (7500 . The results are 



© 0000 RAS, MNRAS 000, 000-000 



Radio loud AGN and the Lx — cr relation of galaxy groups and clusters 7 




100 1000 



a(kms ') CT(kms ') 

Figure 5. The Lx — c relation of radio-loud (left panel) and 
radio-quiot (right panel) subsamples of clusters. The solid line in 
each panel is the best linear fit [equations l[9j and lllOl l]. while the 
dotted line shows the fit to the sample as a whole for comparison 
[equation ||8]l]. 



logL44 = (4.48 ±0.61) log (7500 + (-0.535 ±0.058), (9) 
for radio-loud clusters and 

logL44 = (4.50 ± 0.44) log O-500 ± (-0.530 ± 0.056), (10) 

for radio-quiet clusters. These two relations are plotted as 
solid lines in figure [5] The relation for the cluster sample as 
a whole [Equ. (|8])] is shown as a dotted line in each panel 
for comparison. Again, the fitting parameters are listed in 
Table [1] 

Even though radio-loud clusters are more frequently de- 
tected in the X-ray than radio-quiet objects, the two rela- 
tions in Figure 4 are similar. As we demonstrate in the next 
section, this is simply a selection effect. The majority of indi- 
vidually detected clusters are just above the X-ray detection 
limit in both samples. Since the redshift distributions of the 
two types of cluster are similar (Fig. 1), this forces the mean 
Lx — c" relations of detected objects to be nearly the same. 



4.2 Results for the stacks 

We stack the radio-loud and the radio-quiet clusters inde- 
pendently. Clusters with unclear radio properties are ex- 
cluded from this analysis. 

Radio-loud clusters with a in the range 300 — 900km s~^ 
are stacked in velocity dispersion bins with width 100km s~^ . 
Smaller groups with a < 300km s~^ and bigger clusters with 
a > 900km s~^ are split into two separate velocity dispersion 
bins. Our sample of 134 radio-loud clusters then splits into 
8 stacks. 

As described above, we create a control sample of 134 
radio-quiet clusters selected to have the same a and z dis- 
tributions as the radio-loud sample. This control sample is 
generated from the full sample of 433 radio-quiet clusters by 
picking the radio-quiet cluster that is most similar in z and a 
to each of the radio-loud clusters. These radio-quiet clusters 
are then stacked in exactly the same way as the radio-loud 
clusters. More specifically , for i-th radio-loud cluster that is 
included in the fe-th radio-loud stack, the corresponding i- 
th control radio-quiet cluster is stacked into fc-th radio-quiet 
stack. X-ray detections are obtained for all 16 stacks. 



4-2.1 Surface brightness profiles of the clusters 

We show radial surface brightness profiles for each of our 
eight different velocity dispersion stacks in figure |6] The 
radio-loud clusters are plotted as triangles while the control 
radio-quiet clusters are plotted as squares. Surface bright- 
ness is given in units of photon counts per unit area per sec- 
ond and is plotted as a function of the scaled radius R/Rioo- 
The error in the surface brightness in each radial bin is es- 
timated from the Poisson fluctuations of the photon counts. 
As we can see, after stacking, the S/N of the surface bright- 
ness profiles in all the velocity dispersion bins is sufficiently 
good to enable model fitting to be carried out. 

We use a (3 model [Equ. ^] to fit the surface brightness 
profile for each stack. To reduce the number of free param- 
eters, we fix /3 = 2/3 and estimate Rc and 7o by minimizing 
. We show the best fitting results as solid and dashed lines 
in each panel for radio-loud and radio-quiet clusters respec- 
tively. The best fitting value of Rc for each profile is quoted 
as a label in each panel of Fig. [S] Except for the two lowest 
velocity dispersion bins, radio-loud clusters have more con- 
centrated luminosity profiles (smaller Rc) than radio-quiet 
clusters. 

We show the best fit values of Rc as a function of veloc- 
ity dispersion in figure[7l Radio-loud and radio-quiet clusters 
are represented by triangles and squares, respectively. The 
median Rc of the radio- loud clusters is 0.14J?200, whereas the 
median for radio-quiet clusters is 0.18i?2oo- These two me- 
dian values are shown as horizontal lines in Fig. [T] The core 
radii of our stacked clus ters are consistent wi t h the studies 
of individual clusters bv lNeumann fc ArnaudI (|l999l ). These 
authors found Rc ~ 0.1 — 0.2i?2oo for clusters with /3 ~ 2/3. 
We note that uncertainties in the centroids of individual 
clust er will broaden the X-ray core of the stacked profile 
(e.g. iDai et al.|[2007l ). This effect is not significant in our 
study, however, since the centroids of most of our individual 
detected clusters (which contribute the bulk of the flux of 
the stacks) were identifled before stacking(see Section [3.11) . 
We flnd that if we stack only the clusters with individual de- 
tections, the median Rc values are 0.14i?2oo and 0.17J?200 for 
radio-loud and radio-quiet clusters respectively. The point- 
spread function of the ROSAT telescope also broadens the 
estimated core radii of the clusters, of course, particularly for 
low velocity dispersion clusters which typically have small 
angular size. 



4-2.2 Lx — o relation for the stacks 

In flgureO we show the weighted average X-ray bolometric 
luminosity Lx,s [Equ. (|6])] as a function of the average veloc- 
ity dispersion as for our stacks of radio-loud and radio-quiet 
clusters. We use the BCSE orthogonal regression method to 
fit linear relations between log Lx,s and log erg weighted by 
the errors on both Lx,s and (jg. The error on erg is esti- 
mated from the error on the mean value of o of the stacked 
clusters. Our result is 

log L44 = (4.40 ± 0.53) log (J500 ± (-0.600 ± 0.099) (11) 
for radio-loud clusters and 

log L44 = (4.07 ± 0.24) log (J500 ± (-0.935 ± 0.049) (12) 
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Figure 6. Surface brightness profiles for stacks of clusters in 8 velocity dispersion bins. Triangles show results for radio-loud clusters 
while squares are for the control sample of radio-quiet clusters. The solid and dashed lines show best fit /3 models for radio-loud and 
radio-quiet clusters respectively. The dotted curves show the contribution to the X-ray emission predicted to come from the radio AGN 
itself (see the text in Section |5. II for details). 



for radio-quiet clusters. These relations are shown as solid 
lines in the left and right panels of Fig. |S] for radio-loud 
and radio-quiet clusters respectively. The parameters of the 
fits are also listed in Table [T] The Lx ~ o" relation for the 
individually detected clusters [Equ. (|SJ] is shown as a dotted 
line in each panel for comparison. The Lx — o relation of 
the stacks of radio-loud clusters is very close to that for all 
individual detections (dotted line). This is a coincidence. If 
we stack only the individually detected clusters, we find a 
significantly higher Lx at given cr than predicted by Equ. 
([HI (the dotted line). This is because, at given cr, the mean 
Lx of the stacked objects is higher than their median Lx- 



The latter is what is approximated by a linear fit in log cr - 
log Lx space. 

The slopes of the Lx — o relations for the radio-loud 
and (control) radio-quiet clusters are consistent within the 
1-cr error, but their zero-points differ significantly. At given 
velocity dispersion, the average X-ray luminosity of radio- 
loud clusters is systematically higher than that of radio-quiet 
clusters. To demonstrate the significance of this effect, we fix 
the slope of both relations to be 4.17 [the error-weighted 
mean of the slopes in Equ. (|lip and (I12|l ]. and then re- 
estimate their zero-points. The results are (—0.591 ± 0.031) 
for radio-loud clusters and (—0.937 ± 0.044) for radio-quiet 



© 0000 RAS, MNRAS 000, 000-000 



Radio loud AGN and the Lx — cr relation of galaxy groups and clusters 9 



5 DISCUSSION 



0.4 



0.3 



0.2 



0.1 



200 



400 



600 
a(km/s) 



800 



1000 



Figure 7. The best fit value of Rc for surface brightness profiles 
of clusters stacked in 8 bins of velocity dispersion. Triangles and 
squares represent radio-loud and radio-quiet clusters respectively. 
Rc = 0.18i?200 (solid line) and Rc = 0.14i?200 (dotted line) are 
the median values of Rc for the radio-loud and radio-quiet clusters 
respectively. 
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Figure 8. The average X-ray luminosity Lx g as a function of 
average velocity dispersion as for stacked clusters. The left and 
right panels show results for radio-loud and radio-quiet clusters 
respectively. The solid lines show the best linear fits to the relation 
between log Lx,s and logo-g [Equ. ITTI l and Equ. 1 1121 ]. The Lx — 
a relation for all clusters with individual X-ray detections [Equ. 
||8ll] is shown as a dotted line in each panel for comparison. 



clusters. At given velocity dispersion, the X-ray luminosity 
of radio- loud clusters is on average 2.2 times higher than 
that of radio-quiet clusters. The difference is significant at 
6.4(T and is consistent with our earlier result that the frac- 
tion of radio-loud clusters with individual X-ray detections 
is substantially higher than the corresponding fraction for 
radio-quiet clusters (Section 14.1.11) . The difference is also 
comparable to the difference in normalization b etween cool- 
ing co re and non-cooling core clusters found bv lChen et al.l 
l|2007l ) using the Lx ~ M relation (M is cluster mass) . 



The results presented above demonstrate that the X-ray 
properties of the fCM are correlated with the radio proper- 
ties of the central BCGs. Clusters with radio-loud BCGs are 
more frequently detected in X-ray images. When we stack 
radio-loud and radio-quiet clusters with similar velocity dis- 
persions and redshifts, we find that radio-loud clusters have 
more concentrated surface brightness profiles and higher av- 
erage X-ray luminosities than their radio-quiet counterparts. 

Up to now, we have not considered X-ray emission from 
the central radio AGN itself. If this emission were compara- 
ble to the X-ray emission from the ICM, all our results might 
be explained without any need to invoke a correlation be- 
tween the radio AGN and the state of the intracluster gas. 
We estimate the X-ray emission from the AGN itself in Sec- 
tion EH] 

Furthermore, we have not yet considered how X-ray 
and radio properties correlate with the stellar properties 
of BCGs. As shown by Best et al.(2005, 2007), radio- 
loud AGN occur more frequently in higher stellar mass 
BCGs, which also tend to be found in more X-ray lumi- 
nous clusters (|Edge fc Stewart! Il99ll ). The correlation be- 
tween BCG stellar mass and cluster X-ray luminosity is 
two- fold. On the one hand, more massive clusters typi- 
cally host higher mass BCGs(see top panel of Fig. llUp . On 
the other hand, at given cluster mass (velocity dispersion), 
clusters with higher m ass BGCs tend to be m ore regular 
and more concentrated (|Bautz fc MorganI Il970l), and such 
clust e rs tend to have high er X-ray luminosities (|David et al.l 
1 19991 : iLedlow et all liooa ). This morphology - BCG mass 
depe ndence a ppears related to the presence of the cool 
cores (|Edge fc Stewart, 199L ). Thus, we may expect the mean 
stellar masses of the BCGs to differ between the radio-loud 
and radio-quiet clusters we have stacked above, and we may 
expect clusters stacked as a function of BCG mass to show 
similar differences in X-ray properties to our radio-loud and 
radio-quiet samples. In Section [5.21 we compare the stellar 
properties of the BCGs in our radio-loud and radio-quiet 
samples, and we stack clusters also as a function of BCG 
properties. 



5.1 The X-ray emission from radio AGN 

The X-ray luminosities of radio-loud AGN Lx,agn 
are correlated with t h eir radio luminosit ies La (e.g. 
iBrinkmann et al.1 |2000| : iMerloni et al] |2003| ). Using the 
Lx,AGN — Lr correlation for radio-loud AGN from the study 
of Brinkmann et al. (2000), we estimate X-ray luminosi- 
ties for our sample of radio AGN {Lx,agn) in the ROSAT 
0.1 — 2.4keV band. Here we assume that the distributions 
of radio flux and of X-ray photon energy are power laws , 
•S'm ~ and N{E) ~ E''^ (see IBrinkmann et al-llioool ). 

We then compare our estimated values of Lx,agn with the 
0.1 — 2.4keV X-ray luminosities Lx,icm measured for our 
clusters within the radius Rx- The results are shown in fig- 
ure O The top panel shows results for individually detected 
clusters; the bottom panel shows corresponding results for 
the stacked clusters. For the stacks, we estimate Lx,agn as 
the weighted average of the estimated X-ray luminosities of 
the central radio AGN. 

As can be seen, the estimated X-ray luminosities of the 
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Figure 9. Comparison of the predicted X-ray emission from 
radio AGN{Lx,agn) to the total observed cluster emission 
(Lx icm) for our sample of radio- loud clusters. The AGN lu- 
minosity is estimated from the Lx — Lr relation of Brinkmann 
et al.(2000). The top panel gives results for individually detected 
clusters, the bottom panel for our eight stacks of clusters. The 
dotted lines represent the relation Lx^agn = 0.05Lx,/CA/- 



Figure 10. Median stellar mass (upper panel) and concentration 
(lower panel) for BCGs as function of the velocity dispersion of 
their host cluster. The red triangles and blue squares represent 
radio-loud and control radio-quiet BCGs, respectively. Horizontal 
error-bars show the range of the velocity dispersion bins, while 
vertical error-bars indicate the 32 to 68 percent ranges of the 
distributions. 



radio AGN are typically a small fraction (< 5 percent) of 
the total measured X-ray emission. The enhancement we 
measure in the X-ray luminosity of radio-loud clusters is 
more than a factor of 2 (Section 14. 2[) . Thus, contamination 
of the X-ray emission by the radio AGN is at most a small 
part of the effect we measure. 

One might also ask whether X-ray emission from the ra- 
dio AGN might explain the more concentrated X-ray surface 
brightness profiles seen in figure |6l The dotted curve in each 
panel shows the predicted contribution of the radio AGN 
to the X-ray surface brightness profile. The Lx — Lr rela- 
tion of iBrinkmann et al.1 (|2000l ) was again used to estimate 
the X-ray luminosity of each radio AGN and a Gaussian 
PSF with FWHM = 1 arcmin was used to predict its count 
rate profile. These count rate profiles were then scaled and 
stacked in order to calculate the average AGN contribution 
to the profile of the stack. As we can see, this contribution 
is negligible even in the lowest velocity dispersion bin. 

5.2 The stellar properties of BCGs 

In this section we investigate the relation between our re- 
sults and the stellar properties of our BCGs, as charac- 
terised by their stellar mass and their concentration. The 
concentration is defined as c = Rgo/Rso- where Rqq and 
R50 are the radii including 90 and 50 percent of the fiux 
from a gala xy. It can be used as an indicator of galaxy 
morpholoery (|Shimasaku et al.ll200ll ). For early- type galaxies, 
the stellar mass correlates closely with the mass of the cen- 
tral supermassive black hole. 



We first divide our radio-loud clusters into 8 velocity 
dispersion bins containing equal numbers of objects. The 
median stellar mass and concentration are then calculated 
for both radio-loud and control radio-quiet clusters in these 
8 bins. In figure [10] we show these median values as a func- 
tion of velocity dispersion. For both properties red triangles 
represent radio-loud clusters, while blue squares represent 
the radio-quiet clusters. The horizontal error-bars show the 
range of velocity dispersion for each bin, while the vertical 
error-bars link the 32 and 68 percentiles of the distribution 
in concentration and stellar mass. 

Clearly, the stellar masses of BCGs correlate with the 
velocity dispersions of their host clusters; higher velocity 
dispersion clusters tend to have higher stellar mass BCGs. 
In addition, at given velocity dispersion, the stellar mass 
of radio-loud BCGs is systematically higher than that of 
radio-quiet objects. This mirrors the earlier finding of Best 
et al.(2005, 2007) that higher stellar mass galaxies, both 
BCGs and non-BCGs, are more likely to host radio-loud 
AGN, but that at fixed stellar mass the BCGs are more 
likely to be radio-loud than the non-BCGs. Concentration 
values show no consistent trends with velocity dispersion or 
radio activity and are typically c ~ 3, in the range expected 
for a galaxy with a de Vaucouleurs profile. 

Given that the radio-loud AGN are biased towards 
higher stellar mass galaxies, it is natural to ask whether 
the higher Lx of our radio-loud clusters simply reflects the 
larger stellar masses of their BCGs. To investigate this is- 
sue, we first study the dependence of the Lx — cr relation on 
BCG stellar mass, independent of BCG radio activity. 
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As before, we first divide our 625 clusters of galaxies 
into 8 velocity dispersion bins. We take the range of each 
bin to be the same as in Fig. [8l We then divide the clusters 
in each bin into three equal sub-samples according to the 
stellar mass of their BCGs. Combining velocity dispersion 
bins, gives us three sub-samples, which we refer to as the low, 
intermediate and high mass BCG samples. The median stel- 
lar mass of BCGs in these three samples are 7.76X 10^°A'/q, 
1.35 X 1O"M0 and 2.14 x lO^Af©. For comparison, the 
median stellar mass of the BCGs in our radio-loud and 
control radio-quiet cluster samples are 1.66 X 10" M© and 
1.35 X 10^^ Mq, respectively. Thus, the high BCG mass sam- 
ple has even larger median BCG mass than our radio-loud 
sample, while the low BCG mass sample has even lower me- 
dian BCG mass than our control radio-quiet sample. The 
intermediate BCG mass sample have nearly the same me- 
dian BCG mass as the control radio-quiet sample. 

For each BCG mass sub-sample, we then combine all the 
clusters in each velocity dispersion bin into a single stack. 
For the low BCG mass sample, the X-ray detection of the 
velocity dispersion bin a > 900 kms~^ is not significant, so 
we combine all the clusters with a > 800 kms~^ into one 
stack, thus ending up with seven stacks for the low BCG 
mass sample and eight stacks for the intermediate and high 
BCG mass samples. We show the resulting Lx — o" relations 
in the three panels of Fig. 1111 As before, we use the BCSE 
orthogonal regression method to fit linear relations between 
logLx and logcr weighted by the error on both Lx and a. 
The fitting relations are shown as solid lines each panel of 
Fig. 111! and the corresponding fitting parameters are listed 
in Table 1. The Lx — cr relation for our stacks of radio- loud 
clusters [Equ. (|11|| ] is shown as a dotted line in each panel 
for reference. 

The X-ray luminosity of the clusters shows a very in- 
teresting dependence on the stellar mass of the BCGs. The 
high and intermediate BCG mass samples have have Lx — cr 
relations which are both essentially identical to that of the 
radio-loud clusters. The low BCG mass sample, on the other 
hand, has an Lx — o relation which is significantly different. 
Clusters in the higher velocity dispersion bins are substan- 
tially less luminous than those in the other two samples. 
This results in a much shallower slope, 2.88 ±0.17, than the 
value ~ 4 found in all the other cases. 

The fractions of radio-loud clusters are, as expected, 
higher in the samples with higher mass BCGs. The num- 
bers for the low, intermediate and high mass samples are 
12, 24 and 29 percent, respectively. We note that radio-loud 
fraction in the intermediate and high BCG mass samples do 
not differ very much, while the low BCG mass sample has 
significantly fewer radio-loud objects. Together, all these re- 
sults imply that the stellar mass of BCGs is strongly coupled 
both with their radio activity and with the properties of the 
ICM. We now explore this connection further by examining 
the BCG stellar mass dependence of the Lx — o" relation 
separately for radio-loud and radio-quiet clusters. 

We split our radio-loud sample into two sub-samples 
with similar velocity dispersion and radio properties, but 
with systemically different BCG stellar masses. In detail, 
we separate the 134 radio-loud clusters into 67 pairs with 
very similar velocity dispersion and BCG radio luminosity. 
We then take the clusters with the higher BCG stellar mass 
in each pair to build the high BCG mass sub-sample, and the 




Figure 11. The Lx — c relation for stacks of clusters with high, 
intermediate and low BCG stellar masses (top left, top right and 
bottom left panels respectively). The solid lines in each panel 
show the best linear fits to the relation between log Lx and log a. 
The Lx — o" relation for the stacks of radio-loud clusters [Equ. 
jlTJ] is shown as dotted line in each panel for comparison. 



clusters with lower BCG stellar mass to build the low BCG 
mass sample. We split the 67 high BCG mass clusters into 
6 velocity dispersion bins and make 6 stacks. Clusters with 
a in the range 300 — 600 kms~^ are split into 3 equal bins 
with width 100 kms"^ clusters with 600 < cr < 800 kms~^ 
make up another stack. The remaining groups with a < 
300 kms~^ and clusters with a > 800 kms~^ make up the 
two remaining stacks. For the low BCG mass sample, we 
make 6 corresponding stacks. For each pair of radio-loud 
clusters, if the high BCG mass object is binned into jth 
stack, the low BCG mass object is binned into jth stack 
also. 

We also split the radio-quiet clusters into two sub- 
samples with differing BCG stellar mass. As for the whole 
sample, we first bin the clusters into 8 velocity dispersion 
bins. In each velocity dispersion bin, we then separate the 
clusters into equal low BCG mass and high BCG mass sub- 
samples. For each sub-sample, we combine all the clusters in 
each velocity dispersion bin into a single stack. For the low 
BCG mass sample, we combine the clusters in the two high- 
est velocity dispersion bins (800 kms~^ < cr < 900 km s^^ 
and cr > 900 kms~^), since the significance of the X-ray 
detection of the cr > 900 kms~^ bin is too low. 

We obtain significant X-ray detections for all stacks in 
these four samples. The resulting Lx — cr relations are shown 
in figure [T21 The two sub-samples of radio-loud clusters are 
shown in the top two panels, with the high BCG mass sample 
on the left and the low BCG mass sample on the right. The 
two radio-quiet sub-samples are shown in the same order 
in the bottom two panels. We use the BCSE orthogonal 
regression method to fit linear relations to the logLx and 
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Figure 12. The Lx — t relation for stacks of low BCG mass 
and high BCG mass sub-samples of radio-loud and radio-quiet 
clusters. The results for radio-loud and radio-quiet clusters are 
shown on the top and bottom rows, respectively. High BCG mass 
sub-samples are shown on the left and low BCG mass sub-samples 
on right. The solid lines in each panel show the best linear fits 
to the relations between log Lx and log u. The Lx — c relation 
for stacks of the full sample of radio-loud clusters [Equ. mill ] is 
shown as dotted line in each panel for reference. 

log a data in all these panels, weighting by the errors in both 
Lx and a. The fitting relations are shown as solid lines in 
each panel and the corresponding parameters are listed in 
Table 1. The Lx — <y relation for the radio-loud stacks of 
the full sample [Equ. is again shown as a dotted line 

in each panel for reference. 

Three of these four subsamples show very similar Lx—(y 
relations, all of which are consistent with the relation for 
the full radio-loud sample. Only the radio-quiet subsample 
with low BCG mass shows a significantly different relation 
which in turn is very similar to that shown in Fig. [TT] for 
the subsample with lowest BCG mass, independent of ra- 
dio activity. Thus it seems that the Lx — cr is similar for all 
clusters except those for which are both radio-quiet and have 
an unusually low mass BCG. Among radio-loud clusters we 
detect no dependence of the relation on BCG stellar mass. 
The difference between radio-loud and radio-quiet clusters 
appears to be due entirely to the presence of a subsam- 
ple of relatively high velocity dispersion radio-quiet clusters 
which have both low BCG stellar mass and low X-ray lu- 
minosity. These might plausibly be systems which have not 
yet fully collapsed, so that both their X-ray luminosity and 
their BCG mass are typical of those usually found in lower 
mass relaxed systems. 

To look more carefully for off-sets between our various 
subsamples we have refit the data plotted in Figs [TTI and [T2l 
fixing the slopes a to be 2.97 for the low BCG mass and 
low BCG mass radio-quiet panels and 4.40 for the other five 
panels. The resulting zero-points 6' are listed with their for- 



mal uncertainties in Table 1. The relations in each group 
are formally consistent within their errors, so we see no 
clear evidence justifying a more complex interpretation of 
the data. The most significant and suggestive difference is 
that the full radio-loud cluster sample is about 30% more 
X-ray luminous than the sample of high BCG mass radio- 
quiet clusters, even though the latter sample has a slightly 
larger median BCG mass. This result is significant at just 
over the 2a level, suggesting that even at fixed cluster ve- 
locity dispersion and BCG mass, radio sources prefer to live 
in more X-ray luminous clusters. 



6 CONCLUSIONS 

We have used the ROSAT All Sky Survey to study the X-ray 
properties of a sample of 625 groups and clusters of galaxies 
selected from the Sloan Digital Sky Survey. We focus on 
the Lx — o relation of these groups and clusters, and we 
study whether this relation depends on the radio properties 
of the central galaxy (BCG). A cluster is termed 'radio-loud' 
if its central BCG is a radio-loud AGN, and 'radio-quiet' 
if the BCG is not detected at radio wavelengths. We find 
that the fraction of clusters with individual X-ray detections 
depends strongly on whether the BCG is radio-loud. Radio- 
loud clusters are detected more frequently than radio-quiet 
clusters of the same velocity dispersion and redshift. 

The Lx — o" relations for individually detected radio- 
loud and radio-quiet clusters are very similar, but this is 
purely a selection bias. The majority of detected clusters 
are just above the X-ray detection limit in both samples. 
Since the redshift distribution at each velocity dispersion is 
similar for the two types of cluster, the mean relations for 
detected objects are forced to be similar. 

By stacking the X-ray images of clusters with similar ve- 
locity dispersion, we studied the average X-ray luminosities 
and surface brightness profiles of our clusters as function of 
velocity dispersion. The average X-ray luminosities of radio- 
loud clusters are systematically higher and their luminosity 
profiles are more concentrated than those of radio-quiet sys- 
tems. X-ray emission from the radio AGN itself is by far 
insufficient to explain this doubling of the X-ray luminos- 
ity. Our results demonstrate convincingly and quantitatively 
that the X-ray properties of the intracluster gas correlate 
with the presence of a central radio source in the way first 
suggested by Burns (1990). 

The stellar masses of the BCGs also correlate with their 
radio properties; radio-loud clusters tend to mave more mas- 
sive BCGs than radio-quiet clusters of the same velocity 
dispersion. Those clusters of given velocity dispersion which 
have unusually low mass BCGs tend to be both underlumi- 
nous in X-rays and radio-quiet. This effect is particularly 
pronounced for large velocity dispersion clusters. Among 
radio-loud clusters we find no dependence of X-ray luminos- 
ity on BCG stellar mass. These results can be summarised 
by saying that the only clear dependence of X-ray properties 
on BCG properties is that high velocity dispersion clusters 
in which the BCG is both low mass and radio-quiet tend to 
be several times less X-ray luminous than otherwise similar 
clusters in which the BCG is massive and/or radio- loud. 

Clearly a high central hot gas density is needed for ef- 
fective fuelling of the radio source. If the radio activity is 
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in turn able to heat the surrounding gas and cause its re- 
expansion, the feedback cycle needed to control the growth 
of the central galaxy could be established. The conditions 
which lead to effective AGN fuelling arc clearly related to 
the presence of a massive central galaxy, although the exact 
causal relationship between BCG growth, AGN activity and 
ICM structure remains to be clarified. 
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APPENDIX A: THE CLUSTERS OF GALAXIES INDIVIDUALLY DETECTED IN RASS 

Here we list the basic properties of 142 clusters of galaxies individually detected in RASS. The description of the columns is 
as follows. 

- column 1: the ID of cluster in SDSS C4 cluster catalogue! http://www. ctio. noao.edu~chrism/C47] ) 

- column 2: the Right Ascension(J2000) of the BCG in decimal degrees 

- column 3: the Declination(J2000) of the BCG in decimal degrees 

- column 4: the Right Ascension(J2000) of the X-ray centre in decimal degrees 

- column 5: the Declination(J2000) of the X-ray centre in decimal degrees 

- column 6: the redshift of the BCG 

- column 7; the virial radius -R200 in unit of arcmin 

- column 8: the velocity dispersion in unit of kms~^ 

- column 9: the flux within Rx in the energy band 0.5-2.0kev, in unit of 10^^^ erg s^^cm~^ 

- column 10: the X-ray extension radius Rx, in unit of arcmin 

- column 11: the fractional error on count rate(flux) within Rx 

- column 12: the extension correction factor /e [equation (|4}] 

- column 13: the bolometric X-ray luminosity after the extension, in unit of 10** erg s~^ 

- column 14: the classification of radio properties of BCG, 1 for the clusters with BCG to be radio-loud AGN, for the 
clusters with BCG to be radio-quiet, —1 for others 
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Figure Bl. The RASS image of tlie C4 cluster 1094. Tlie cross sliows tiie position of the BCG and the big circle enclosed the region of 
the cluster within a distance -R200 (~ 8 arcimin) from the cluster centre. The small circle shows the X-ray point source detected near the 
cluster centre. The angular distance from the X-ray source to the position of the BCG is ~3 arcmin. This X-ray source is judged to be 
a contaminant and not associated with the X-ray emission of the cluster. 




Figure B2. The RASS image of the C4 cluster 1049. The cross shows the position of the BCG and the small circles show the ML 
detected X-ray sources. The inner big circle shows the region within a distance -R200 (~ 16.5 arcimin) from the cluster centre, while the 
outer circle has a radius of ii200 + 6 arcmin. The area inside the ring from -R200 to R200 + 6 arcmin is used for the background estimation. 
In this image, sources 48 and 55 are considered as contaminants, while sources 52 and 53 are considered as part of the extended X-ray 
emission from the cluster. 
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APPENDIX B: EXAMPLES OF OUR PROCEDURE FOR DETECTING CLUSTERS IN THE RASS 

In the following, we show two example images where decisions about the X-ray detections had to be made by eye. The first 
example deals with the process of determining of the X-ray centre of the cluster (Section [3TT}. Some X-ray sources, which are 
although inside the 5 arcmin aperture of the BCGs, are judged as contaminating sources rather than the X-ray centres of 
the clusters (Fig. Bl). The second example deals with the procedure of distinguishing contaminating X-ray sources from the 
extended X-ray emission of the cluster (Fig. B2). 
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